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1.  Introduction 

Quantum  communication  is  a  laser  communication  technology  that,  in  addition  to  very 
high  data  rate  and  low  power  requirements  of  the  transmitters,  offers  unprecedented  data 
security.  Optical  communication  in  general  is  very  popular  when  high  security  is  important 
because  inherently  small  beam  divergence  angles  facilitate  low  probability  of  interception  and 
low  probability  of  detection  (LPI/LPD).  However,  when  additional  immunity  to  eavesdropping  is 
required,  data  encryption  may  be  necessary. 

Optical  communication  offers  the  unique  feature  of  quantum-based  encryption  due  to  the 
inherent  properties  of  light  used  as  a  carrier  signal.  Current  research  efforts  are  aimed  at  using 
various  quantum  states  to  perform  data  encoding;  however,  polarization-based  techniques  are 
still  the  most  popular  ones  for  a  variety  of  tasks,  including  quantum  communication  (QC), 
quantum  key  distribution  (QKD),  and  keyed  communication  in  quantum  noise  (KCQ). 

For  many  practical  needs,  quantum  communication  systems  must  support  operation 
between  mobile  platforms,  which  hinges  upon  several  innovations.  In  particular,  successful 
pointing,  acquisition,  and  tracking  (PAT)  require  the  use  of  a  beacon  signal  and  the  capability  of 
accurate  and  agile  alignment  of  the  line-of-sight  (LOS)  between  the  communicating  terminals 
performed  over  a  large  field  of  regard.  While  mechanical  devices,  such  as  gimbals,  offer 
relatively  slow  tracking  over  a  very  wide  range,  they  lack  in  pointing  bandwidth  necessary  for 
rejecting  high  frequency  vibrations  and  beam  deflection  caused  by  the  optical  turbulence.  In 
contrast,  fast  steering  and  especially  non-mechanical  devices,  such  as  Bragg  cells,  enjoy  very 
high  bandwidth  (on  the  order  of  several  kHz),  but  their  effective  range  is  very  small.  Inherent 
limitations  of  both  gimbals  and  fast  steerers  result  in  shortcomings  of  the  entire  PAT  system 
when  either  of  these  devices  is  used  as  a  sole  beam  steerer.  Therefore,  focus  needs  to  be  shifted 
to  hybrid  architectures,  exploiting  the  advantages  of  the  constituting  elements. 

The  proposed  research  will  utilize  a  mechanical  gimbal  (such  as  Omni-Wrist  or  another 
commercially  available  device)  for  extended  range  of  optical  connectivity,  and  a  fast  beam 
deflector  to  create  a  hybrid  beam  steering  system  capable  of  exercising  a  very  high  positioning 
bandwidth  over  a  full  hemisphere  of  steering  angles.  System  design  process  will  include  the 
solution  of  such  underlying  problems  as  the  development  of  the  mechanical  and  optical 
subsystems,  mathematical  description  of  the  hybrid  device,  optimal  task  distribution  between  the 
mechanical  and  non-mechanical  positioning  components,  and  coordination  of  the  operation  of 
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the  “coarse”  and  “fine”  system  controllers.  The  efficiency  of  the  developed  system  under  various 
operational  conditions  will  be  investigated  and  compared  against  known  designs.  It  is  proposed 
to  develop  advanced  control  strategies  that  would  assure  a  highly  coordinated  operation  of  both 
system  components,  thus  resulting  in  a  beam  steerer  with  previously  unknown  range,  agility  and 
accuracy. 
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2.  Gimbal  Device  for  Wide-Range  (Coarse)  Beam  Steering 
Omni-Wrist  III  (see  Fig.  2.1)  is  proposed  as  a  possible  device  for  coarse  steering.  It  is  a 
new  sensor  mount  developed  under  Air  Force  funding  that  emulates  the  kinematics  of  a  human 
wrist.  Driven  by  two  linear  motors  and  computer  controlled,  it  is  capable  of  a  full  180° 
hemisphere  of  pitch/yaw  motion.  A  comprehensive  laboratory  testing  of  one  of  few  existing 
devices  of  this  type,  installed  in  the  Laser  Communications  Research  Laboratory  at  Binghamton 
University,  has  resulted  in  the  establishment  of  a  complete  transfer  matrix-type  model  relating 
pitch/yaw  coordinates  of  the  sensor  mount  to  the  motor  encoder  signals. 


Figure  2.1.  Omni-Wrist  III  Sensor  Mount 

In  contrast  to  traditional  robotic  manipulators,  the  actuators  driving  Omni-Wrist  III  are 
not  located  in  the  joints  but  rather  attached  to  the  links,  mimicking  the  attachment  of  muscles  to 
bones  in  biological  structures.  The  resultant  device  is  a  two-degree-of-freedom  system  capable  of 
a  full  180°  hemisphere  of  singularity-free  yaw/pitch  motion  with  up  to  5  lbs  of  payload.  In 
comparison  to  traditional  gimbals  positioning  devices,  Omni-Wrist  III  enjoys  increased 
bandwidth  due  to  a  greater  power/mass  ratio,  and  reduced  inertia  and  friction.  However,  its 
mechanical  design  does  not  eliminate  nonlinearities  and  cross-coupling,  complicating  the 
controls  task.  In  our  previous  work  we  developed  the  solution  to  the  inverse  and  forward  pose 
kinematics  problem,  and  investigated  the  dynamics  of  the  system,  as  outlined  in  the  following 
sections  of  this  chapter.  The  resulting  mathematical  model  of  Omni-Wrist  III  builds  a 
framework  for  the  synthesis  of  advanced  control  strategies  required  to  utilize  this  device  to  its 
full  potential. 
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2.1.  Pose  kinematics. 

The  cross-shaped  moving  sensor  mount  is  connected  to  the  stationary  platform  of  the 
same  shape  through  four  identical  legs,  each  comprising  three  links  and  four  revolute  joints.  Two 
of  the  legs  are  redundant  for  the  development  of  the  pose  kinematics  and  are  omitted  in  the 
kinematic  diagram  (see  Fig.  2.2).  The  position  and  orientation  of  the  joints  is  symbolized  by 
short  thick  lines.  The  linear  motors  are  connected  to  the  two  bottom  links  as  shown  in  Fig.  2.3(a). 


Figure  2.2.  Omni-Wrist  III  Kinematic  Diagram 


bcosQ-d  d 


(a) 


(b) 


Figure  2.3.  (a)  Connection  of  the  Actuator;  (b)  Azimuth,  Declination,  Yaw  (a).  Pitch  (fi) 
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In  the  three-step  process  of  finding  the  pose  kinematics  solutions,  the  correspondence  between 
the  azimuth  and  declination  and  yaw  and  pitch  coordinates  is  found,  as  well  as  between  angles  0\ 
through  Os  and  the  yaw  and  pitch  angles  and  between  the  angles  6\  and  O5  describing  the  rotation 
of  the  joints  connecting  Leg  A  and  Leg  B  to  the  stationary  platform  and  the  position  of  the 
actuators.  The  structure  of  Omni-Wrist  III  is  captured  in  the  kinematic  diagram  in  Fig.  2.2 
showing  the  x  and  z  axes  of  all  intermediate  coordinate  frames  defined  according  to  [1].  The 
corresponding  Denavit-Hartenberg  parameters  could  be  found  to  derive  the  transformation 
between  the  stationary  frame  and  the  sensor  mount  frame  through  Leg  A  and  Leg  B  [2],  Both 
transfonnations  are  equal  to  the  transformation  into  the  yaw,  pitch  and  roll  coordinate  system: 

A1xA2xA3xA4-B0xBlxB2xB3xB4-C,  (2-1) 

where  C  is  defined  as 


c  c 

css  -sc 

CSC  +s  s 

X 

a  0 

saca 

a  0  n  an 

SaSaSn  +  CaCn 

a  0  n  an 

SaSQCn-CaSn 

y 

-s„ 

cnsn 

cc„ 

z 

(2.2) 

0 

0 

0  n 

0 

0  n 

0 

1 

where  indexes  n,  o  and  a  correspond  to  the  roll,  yaw  and  pitch  coordinates,  and  S  and  C  denote  a 
sine  and  a  cosine,  respectively. 

2.1.1.  Inverse  Pose  Kinematics.  For  the  solution  to  the  inverse  pose  problem,  the  actuator 
encoder  values  need  to  be  found,  which  correspond  to  the  given  azimuth  and  declination 
coordinates.  This  is  achieved  in  three  steps.  First,  the  yaw  and  pitch  coordinates  are  found,  which 
correspond  to  the  azimuth  and  declination  coordinates,  followed  by  determining  the  values  of  the 
eight  joint  variables  (angles  0\  through  Os)  corresponding  to  the  yaw  and  pitch  coordinates.  In  the 


last  step,  the  correspondence  between  the  actuator  encoder  values  and  the  values  of  0\  and  O5 
(parameters  of  joints  between  the  bottom  links  and  the  stationary  platform)  is  established. 

The  solution  to  the  first  step  can  be  found  by  investigating  Fig.  2.3  (b),  which  shows  the 
graphical  representation  of  these  equations: 

1  +  tan 2  ( az )  cos 2  (dec) 


cos 


:  (yaw)  ■ 
cos  (pitch) - 


l  +  tan2(az) 
cos(<iec) 


cos 


(yaw) ' 


(2.3) 


(2.4) 
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Application  of  the  second  step  is  explained  in  detail  in  [2],  and  it  results  in  the  following 
relationships  between  two  pairs  of  variables  in  each  leg: 

9 =  9,  +  n 


for  leg  A  and 

9b-91Jrl l 
9&  —  —9S 


(2.6) 


for  leg  B,  simplifying  significantly  the  analytic  solution  for  the  remaining  variables. 

The  key  to  the  last  step  of  the  solution  to  the  inverse  kinematic  problem,  the  connection 
of  the  actuators,  lies  in  Fig.  2.3(a),  which  shows  a  geometric  model  of  the  situation.  The 
corresponding  equations,  which  provide  values  for  Ax\  and  Ax2,  the  actuator  encoder  positions, 
from  6\  and  6*5  are 


(a+b  sin  0X  )2  +  ( b  cos  6X  -  d)2  =  (c  +  Axx  )2 , 
{a  +  b  sin  05  )2  +  ( b  cos  05  -  d)2  =  {c  +  Ax2  )2 , 


where  c  represents  the  extension  of  the  actuator  corresponding  to  encoder  value  of  zero.  The 
model  implemented  in  the  computer  controller  supplied  with  the  Omni-Wrist  III  but  without 
documentation  was  used  as  a  black  box  to  provide  data  for  the  detennination  of  a,  b,  c,  and  d 
utilizing  genetic  optimization  [2],  [3], 

2.1.2.  Forward  pose  kinematics.  Similarly  to  the  inverse  kinematics  case,  the  forward 
pose  problem  is  solved  in  three  steps,  now  in  reversed  order.  At  first,  angles  6\  and  O5  are 
detennined  from  the  knowledge  of  the  actuator  encoder  values  Ax  1  and  Ax2.  Then,  the  remaining 
angles  02  through  6*4,  6()  through  9%  and  the  yaw  and  pitch  are  derived  from  angles  0\  and  02 
followed  by  the  conversion  of  the  yaw/pitch  coordinates  into  azimuth  and  declination.  In  the  first 
step,  the  quadratic  terms  in  (2.7)  are  expanded  to 

a 2  +b 2  +  d2  -  (c  +  Ax,  Y  +  lab  sin(#. )  =  2bd  cos (6, ), 

12  81 

a 2  +b 2  +d 2  -  (c  +  Ax 2  )2  +  lab  sin(6)5 )  =  Ibd  cos(#5 ). 


Squaring  (2.8)  and  rearranging  produces 

sin2  9x4b2[a2  +  d^+ sm9xAab{cr  +b 2  +d 2  ~{c  +  Axl)1^)+{cr  +b 2  +d 2 -(c  +  Axxff  -4b2d2  =0, 
sin2  954b2(a2  +d2')+sin954ab(ci2  +b 2  +d 2  -(c  +  ^4x0)2)+(n2  +b 2  +d 2  -(c  +  Ax^f  )  -4b2 d1  =0, 

which  are  simple  quadratic  equations  and  can  be  easily  solved  to  give  6\  and  <95. 
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In  the  next  step,  S3  is  expressed  from  element  of  the  transformation  matrices  in  [2]  as 


c,  s.  c, 

and  manipulated  further  to  produce  [2] 

c,  ^  +  ^(s,  -  s,  +  ^0,5,  =  Ct . 

L1  ^ a  C1 

Squaring  and  rearranging  (2.1 1)  forms 


(2.10) 


(2.11) 


Cl 


/ 

f  sl  ) 

f 

1+1 

—c5 

+  2  C6 

c, 

V 

v^i  y 

) 

V 

C,  +^-(St  -S5) 


1  -c„ 


f  c  \ 

—C. 

C, 


=  0 ,  (2.12) 


which  can  easily  be  solved  to  obtain  0b.  In  an  equivalent  procedure,  Se  can  be  expressed  from 
element  of  the  transfonnation  matrices  to  produce  [2] 

(2.13) 


C5  ^ {S5  -  S \  +  ^CtS2=-C2, 


' 5  w a  ^5 

which  can  be  squared  and  rearranged  to  form 


Ci 


s  $  ) 

2  / 

1  +  — c. 

+  2  C2 

c 

\  ^5  J 

V 

Cs+^iSs-S,) 


1  -c„ 


( 5  A 

—  Q 
vQ  y 


Q-f-ffo-S,) 


1-C„ 


=  0  (2.14) 


to  obtain  6%. 

In  a  similar  fashion,  the  yaw  and  pitch  coordinates  60  and  0a  can  be  obtained  as  follows 

S,C4(C2C3  -  5253Cj-(C,53C4  ~StS2S4)Sa -ClS4Ca  =  SaCa 

C4{S2C2+C2S2Ca)-C2S4Sa=-S0 

and  converted  into  azimuth  and  declination  by  using  the  following  equations,  which  are 
represented  graphically  in  Fig.  2.3(b) 
cos  (dec)  =  cos(va\i’)cos(pitch) 
tan(crz)  = 


sin  (yaw) 


cos 


(yaw)sin(pitcli) 


(2.15) 


(2.16) 


2.2.  Dynamics. 

Due  to  the  necessity  for  a  successful  development  of  a  control  system,  the  dynamic 
properties  of  the  Omni-Wrist  were  also  investigated.  Because  of  the  complicated  structure  of  the 
device,  a  response  was  collected  at  thirteen  different  points  on  the  hemisphere  covered  by  Omni- 
Wrist  III,  with  azimuth  and  declination  values  of  (0,0),  (0,30),  (60,30),  (120,30),  (180,30), 
(240,30),  (300,30),  (0,60),  (60,60),  (120,60),  (180,60),  (240,60),  and  (300,60).  The  actuators 
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exhibited  the  velocity  response  of  a  first  order  system;  therefore  a  least  squares  estimation 
procedure  was  implemented  to  fit  the  response  of  a  first  order  system  to  the  response  of  the 
Omni-Wrist  in  the  form 


s  +  a 


(2.17) 


The  results  of  the  estimation  are  summarized  in  [2]  and  [3],  where  the  placement  of  poles  at  the 
thirteen  different  locations  for  axis  1  and  axis  2  in  positive  and  negative  direction  were  found. 
All  poles  are  located  in  the  95%  confidence  interval  around  47  for  the  positive  direction  and  25 
for  the  negative  direction  for  both  axes. 

While  the  magnitude  of  the  response  differs  in  the  positive  and  negative  direction,  it  is 
possible  to  model  the  dynamics  of  the  system  using  only  two  transfer  functions  for  each  actuator 
with  voltage  input  and  position  output  as 


G?{s)  = 


2.3-106  2.4  •  106  1.7  •  106  1.9-106 

Gi  (©  =  2  »  G2  W  =  ^ — — ,  G2[s)=- 


sz  +  25s 


sz+47s 


s~  +25s 


(2.18) 


sz+47s 

2.3.  Overall  Model. 

The  developed  mathematical  model  is  a  crucial  starting  point  for  the  design  of  an 
efficient  control  system.  Fig.  2.4  represents  the  proposed  configuration  of  the  mathematical 
model  of  Omni-Wrist  that  comprises  two  modules: 

DYNAMICS  represents  the  dynamics  of  two  independently  operating  linear  actuators  coupled  to 
the  sensor  mount  through  a  series  of  links  and  joints.  It  includes  two  transfer  functions,  Gi(s)  and 
G2(s),  describing  the  typical  linear  relationships  between  the  control  efforts,  voltages  U/  and  l A, 
and  the  resultant  linear  displacements,  x  and  y. 

DYNAMICS  KINEMATICS 


Ui 

U2 


Figure  2.4.  Configuration  of  the  Omni-Wrist  model 

KINEMATICS  describes  the  nonlinear  static  relationship  between  the  linear  actuator 
displacements,  x  and  y,  and  the  angular  displacements  of  the  platform,  0i  and  02.  Elements 
( Pjj(x,y),  i,j=l,2,  reflect  the  complex  kinematics  of  the  Omni-Wrist  structure.  Analysis  of  the 
system  has  resulted  in  a  system  of  trigonometric  equations;  however,  these  equations  are  too 
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complex  for  any  direct  use,  and  it  seems  to  be  more  practical  to  represent  both  the  kinematics 
and  inverse  kinematics  of  the  device  by  a  sequence  of  three  transformations  [2],  Development  of 
this  module  implies  the  solution  of  the  direct  pose  kinematics  problem  utilizing  the  Denavit- 
Hartenberg  approach  and  finding  transformation  of  the  linear  encoder  readings  into  joint 
coordinates,  then  into  the  yaw,  pitch  and  roll  angles,  and  finally  into  the  azimuth  and  declination 
coordinates. 
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3.  Acousto-Optic  Device  for  Agile  (Fine)  Beam  Steering 

In  order  to  transition  the  technology  being  developed  under  various  quantum 
communication  R&D  initiatives  to  military  and  commercial  end  users,  several  outstanding 
problems  must  be  solved.  Among  these  is  the  development  of  agile,  compact  beam  steering.  This 
is  a  requirement  for  Air  Force  applications  where  optical  communications  gear  must  be  fitted  to 
an  airframe,  which  places  severe  demands  on  the  pointing  and  tracking  element. 

The  state  of  the  art  in  beam  steering  employs  gimbals  for  coarse  steering  augmented  by 
fast  steering  mirrors  (FSM)  for  fine  steering.  However,  the  current-generation  mirrors  are  barely 
adequate  for  airborne  applications.  In  spite  of  their  simplicity  and  well  established  usage  of  tip- 
tilt  FSM  as  a  fine  steering  device,  they  have  a  few  well  known  drawbacks:  relatively  high  total 
weight,  their  actuators  require  fairly  large  currents  and  consequently  high  plug-wall  power  of  the 
electrical  drivers.  Inertia  and  mirror's  eigenmode  interaction  typically  limit  their  bandwidth  to 
less  than  1  kHz,  while  it  is  commonly  known  in  atmospheric  laser  beam  propagation  that 
adaptive  systems  faster  than  1  kHz  are  highly  desirable  to  overcome  optical  turbulence  induced 
by  the  airframe’s  skin  and  bow  shock  flows.  Development  of  a  practical  system  demands  a 
replacement  for  the  FSM  that  is  small,  light,  low  power,  and  capable  of  improving  the  control 
bandwidth  by  a  factor  of  4  or  more. 

Acousto-optics  deflection  (AOD)  technology  can  offer  much  larger  bandwidth  (more 
than  20  kHz)  within  roughly  the  same  excursion  range.  We  believe  that  a  system  can  be 
developed  that  is  substantially  smaller  and  lighter  than  the  FSMs,  that  requires  almost  an  order  of 
magnitude  less  power,  and  that  can  achieve  control  bandwidths  of  several  kHz.  As  a  practical 
example,  latest  models  of  handheld  barcode  scanners  are  increasingly  employing  miniature  AOD 
elements,  particularly  for  scanning  of  fast  moving  objects. 

3.1.  Acousto-Optic  Deflection. 

An  acousto-optic  cell  utilizes  the  effect  of  Bragg  diffraction  of  the  laser  beam  incident 
upon  a  volume  grating  (see  Fig.  3.1).  An  ultrasonic  wave  is  used  to  create  regions  of  expansion 
and  compression  inside  the  Bragg  cell,  causing  changes  in  density.  The  index  of  refraction  is 
then  periodically  modulated  and  the  medium  becomes  equivalent  to  a  moving  phase  grating  [4] 

An(z,t)  =  An  sin(wst  -  ksz),  (3.1) 

where:  z  is  the  position  inside  the  Bragg  cell  along  the  vertical  axis;  ws  and  ks  =  acoustic 
frequency  and  wave  number,  respectively. 
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©D  =  A£/nv-2©B 


|  acoustic 
1  transducer 


frequency  (y 


^  =  laser  wavelength 
f  c  =  center  acoustic  freq 

n  =  refractive  index 
v  =  acoustic  velocity 


Figure  3.1.  Bragg  cell  operation 

The  angle  of  incidence  is  selected  in  such  fashion  that  the  conservation  of  energy  and  the 
principle  of  momentum  conservation  between  the  acoustic  and  optical  wave  vectors  during  light- 
sound  interaction  is  preserved  [5].  It  leads  to  a  mathematical  expression  commonly  known  as  the 
Bragg  angle 

0B  =  sin"1  (2fc  / 2nv),  (3.2) 

When  the  acoustic  frequency  applied  to  the  Bragg  cell  is  varied  from  fc  to  (fc+Afs),  there 
is  a  change  in  the  magnitude  of  the  sound  vector  equal  to  AK=2 7i(Afs)/v.  As  a  result,  the 
diffracted  beam  will  propagate  along  the  direction  that  least  violates  the  momentum  conservation 
principle.  This  change  in  the  sound  vector  results  in  a  small  angular  motion  of  the  deflected 
beam  and  is  found  to  be  proportional  to  the  frequency  of  the  input  acoustic  signal  via 

A©  =  (A*  Afs )  !{n  *  v) ,  (3.3) 

Hence  the  direction  of  the  diffracted  beam  could  be  controlled  by  the  frequency  of  the 
acoustic  wave/  with  a  deflection  angle 

0D=JU(nv)*(f-fe)  (3.4) 

A  closed-loop  configuration  of  an  acousto-optic  system  with  two  Bragg  cells  is  shown  in  Fig. 3.2. 
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Vertical  cell  Quac|rant 


3.2.Dynamics  of  Acousto-Optic  Steering. 

The  deflection  model  (3.2)  suggested  as  a  result  of  theoretical  study  of  the  acousto-optic 
phenomena  uses  the  acoustic  frequency  of  the  transducer  as  an  input  and  deflection  angle  of  the 
laser  beam  as  an  output;  however,  this  equation  does  not  reflect  any  transient  when  steering  is 
performed.  Bragg  cells  are  characterized  by  very  fast  beam  steering,  and  the  following  dynamic 
changes  describe  the  process.  The  deflection  angle  Od  changes  while  the  acoustic  wave  traverses 
the  laser  beam  that  has  width  w.  Therefore,  response  time  is  determined  by  the  ratio  of  the  Bragg 
cell  aperture  to  the  acoustic  velocity  in  the  interaction  media  A  typical  value  for  tellurium 
dioxide  lies  in  the  10- lOOps  range;  however,  advanced  controls  are  still  required  for  proper 
operation  of  the  beam  steering  system,  especially  in  the  presence  of  platform  vibrations  or  strong 
atmospheric  turbulence.  Since  the  deflection  angle  does  not  experience  any  overshoot  [6], 
dynamics  of  this  process  can  be  best  described  by  a  first  order  system  (lag  filter).  As  a  result,  the 
model  of  the  Bragg  cell  can  be  established  in  the  form  of  the  following  transfer  function  [7] 

G(s)  =  0D/(f-fc)  =  [A/(nv)]  *  |  wb/(s+  wb)\  (3.5) 

where  wb  is  a  parameter  of  the  lag  filter  modeling  access  time  of  the  Bragg  cell. 

Definition  of  the  parameters  of  the  above  model  and  extensive  study  of  the  acousto-optic 
phenomena  suggests  no  cross-coupling  between  two  sequentially  mounted  Bragg  cells,  when 
two-dimensional  steering  is  required. 
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4.  Hybrid  Beam  Steering  System 

As  outlined  in  the  previous  chapters,  state-of-the-art  beam  steering  devices  cannot  offer 
adequate  combination  of  features  to  satisfy  the  range  and  speed  requirements  imposed  by 
application  of  lasers  in  mobile  communication  systems.  Hence,  the  pointing,  acquisition,  and 
tracking  (PAT)  task  of  any  laser-based,  and  particularly,  quantum  communication  system  is 
superficially  decomposed  into  the  “slow-high-magnitude  motion”  problem,  known  as  coarse 
steering,  and  the  “fast-low-magnitude  motion”,  known  as  fast  steering.  This  research  is  focused 
on  the  development  of  a  hybrid  system  addressing  the  needs  of  the  entire  PAT  task,  utilizing  two 
different  advanced  beam  steering  technologies  fused  together  by  hierarchical  control.  The  first 
technology  is  the  Omni-Wrist  III  mechanical  system  known  for  its  superior  dynamics,  in 
comparison  with  traditional  gimbals,  and  a  full  hemisphere  steering  range  [2].  The  second 
technology  is  the  acousto-optic  Bragg  cell,  virtually  inertia-less  but  with  a  very  limited  steering 
range  [8],  [9]. 

4.1.  Proposed  Approach. 

Fig.  4.1  illustrates  graphically  the  idea  of 
wide  range  connectivity  facilitated  by  a 
mechanical  device  (such  as  the  Omni-Wrist 
gimbal)  combined  with  high  bandwidth  of  a 
narrow-range  agile  steerer  (Bragg  cell).  As  a 
configuration  example,  two  Bragg  cells,  required 
to  provide  two-coordinate  beam  deflection,  could 
be  integrated  into  an  optical  setup  and  placed 


Laser  beam  position 


Figure  4.1.  Range-bandwidth  of  a  hybrid  device 
directly  on  the  sensor  mount  of  the  Omni- 

Wrist,  as  illustrated  in  Fig.  4.2.  One  of  the 
most  attractive  features  of  the  proposed 
system  is  its  ability  to  effectively  cover  a 
wide  area  in  the  “range  -  bandwidth” 
domain  in  the  sense  that  the  high  bandwidth 
capability  could  be  “delivered”  to  any 
location  on  the  hemisphere  of  system 


Figure  4.2.  The  hybrid  steerer  concept 
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operation  (it  should  be  realized  that  high  frequency  components  of  any  steering  task  have  small 
amplitudes). 

4.2.  Omni-Wrist  III  Control  System. 

4.2.1.  Control  Synthesis.  When  used  as  a  steering  device  for  pointing,  acquisition,  and 
tracking  (PAT),  Omni-Wrist  III  could  be  viewed  as  a  two-input-two-output  system  that  positions 
the  laser  beam  over  a  wide  range  of  azimuth  and  declination  angles.  Its  model,  with  the  structure 
shown  in  Fig.2.4,  could  be  presented  with  a  single  nonlinear  transformation  that  accounts  for 
both  kinematic  and  dynamic  properties  as  follows. 


©r 

©a  Ay) 

®i2  Ay) 

G,A) 

0 

gu(s’x>  y) 

g  12  A  x,y) 

y 

(4.1) 
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0 

A  A)_ 
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Hence,  the  output  of  either  dynamic  channel  (azimuth  or  declination)  can  be  found  as 

0,  =  gllU,  +gijUj  (4.2) 

Let  di  (disturbance  signal)  represent  cross-coupling  effects,  nonlinearities  and  other  unmodeled 


dynamics  [10],  such  that  (4.2)  becomes 

=  GiUi  +  dt  =  — 1 LUi  +  di  =  — - Ui  +  dt , 

A  s  ~  +  as 

where  Nt  and  A  -  numerator  and  denominator  polynomials  of  G'„  respectively. 
Then  suggested  control  signal  is  formed  as  follows 

Tt=Ntut  =Dt®,-Dtdt 


(4.3) 


(4.4) 


For  simplicity  the  subscript  i  identifying  the  dynamic  channel  can  be  omitted  resulting  in 
the  following  differential  equation  when  (4.4)  is  presented  in  the  time  domain 

T  =  ®"+a@'+d  ,  (4.5) 

where  a  new  disturbance  term  is  d  =  Didi . 


The  proposed  control  system  for  each  dynamic  channel  has  a  unity  gain  feedback  and 
three  modules:  conventional  controller,  adaptive  feedback  and  feedforward  controllers  in  the 
following  fonn 

T  =  [l0jedt  +  lle\  +  [k0e  +  kle']  +  [q0@r  +  ql@r'+q2®r"],  (4.6) 

where  /,,  and  qf  -  controller  gains. 

System  configuration  for  one  channel  is  presented  in  Fig.  4.3. 
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Figure  4.3.  Decentralized  adaptive  control  system 

Substituting  the  first  term  of  the  above  equation  with/,  as  shown  in  Fig.  3,  and  combining  (4.5) 
and  (4.6)  results  in 

e"+(a  +  kx)e'+kQe  =  d  -  f  - q0Q r  ~(ql  -a1)&r'-(q2  -1)0,."  (4.7) 


The  above  second-order  differential  equation  has  a  matrix-vector  equivalent  that  could  be 
obtained  by  introducing  X=[e  e  ]T,  as  follows 
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(4.8) 


If  vector  Xm  =[em  em]T  represents  the  desired  error  signals,  then  its  dynamics  of 
convergence  towards  zero  can  be  described  by  the  following  state  equation  written  in  the 
canonical-controllable  form. 


X,  „  =  DX„  = 


0 
-  a 


o 


1 

-  a 


X. 


(4.9) 
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Since  (4.9)  is  a  stable  system,  there  exists  a  solution  to  the  Lyapunov  equation  [10] 

PD  +  DtP  =  -Q, 

where  P  and  Q  -  positive  definite  matrices. 

Introduction  of  vector  E=Xm-X  leads  to  a  formal  mathematical  definition  of  the  process  of 
error  convergence 


0  1 

0  1 

E  = 

-  a„  -  a”‘ 

E  + 

k0  -  a™  a  +  kx-  a”' 
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qx-  ci 
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(4.10) 


A  control  law  based  on  a  Lyapunov  function  obtained  from  (4.10)  would  ensure  that, 
given  any  initial  condition,  E  converges  asymptotically,  and;  therefore,  the  actual  error  trajectory 
X  will  track  the  desired  error  trajectory  Xm  that  converges  to  zero.  Consider  the  following 
positive-definite  function  as  a  Lyapunov  candidate. 

V  =  ET PE  +  Q0(f  ~  d)1  +  0,  (k0  -  a ;  )2  +  02 (a  +  *,  -  < )2  + 

+  Qrf02  +QM  i  -«)2  +(?5(92  -!)» 


where  Q,  -  positive  scalars. 

Its  derivative  must  be  negative-definite  in  order  to  claim  that  (4.1 1)  is  a  Lyapunov  function.  For 
obtaining  an  analytical  expression  the  assumption  that  the  controlled  plant  is  “slowly  time- 
varying”  compared  to  the  control  effort  is  suggested  [10];  therefore,  d  =  0,  and  differentiating 
(4.11)  along  the  trajectory  defined  by  (4.10)  results  in 

v  =  -EtQE  +  2  Q0  (/  -  d)f  -if-  d)r  +  2  Qx  (k0  -  a”  )k0  -  (k0  -  )re  + 

2 Q2  (a  +  kl  -  a'”  )k  -  (a  +  kx  -  a”  -  q0r®r  +  (4- 12) 

204 (9i -a)4i  -  (9i  ~a)r®r'+2Q5(q2  -  \)q2  -  (q2  -  \)r®r", 

where 


r  =  w0e  +  wle',  (4.13) 

and  wo  and  wi  are  positive  weighting  coefficients.  Grouping  tenns  of  the  equation  (4.12)  results 
in 

V  =  -EtQE  +  (/  -  d)(2Q0f  -  r)  +  (k()  -  a™)(2Qxk0  -  re)  +  (a  +  kt-  a[n)(2Q2k  -  re')  + ...  ^  ^ 
+  9o(2039o  -'•©,  )  +  (9i  -a)(204?i  ~r®r')  +  (q2  -1)(2 Q5q2 -r@,  ") 

While  there  could  be  multiple  solutions  to  the  control  synthesis  problem  that  result  in  (4.14) 
being  negative-definite,  the  most  natural  way  to  select  the  adaptation  law  is  as  follows. 
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2Qj-r  =  0, 

2Q\k0  -  re  =  0 , 

2Q2kl  -  re'=  0 , 

2Q3q0-r®r=0,  (4.15) 

2Q4qt  -r@r'=  0, 

2Q5q2  -/•©,."=  0, 

Hence,  the  time  derivative  of  function  V  becomes 

V  —  —EtQE.  (4.16) 

Solving  (4.15)  for  unknown  variables  provides  expressions  for  the  conventional  controller 

/  =  d^rdt  =  Sw0^edt  +  Sw{e  =  l0  ^edt  +  lxe  (4.17) 

and  equations  for  the  adjustable  gains  of  adaptive  controllers 
k0  =  ax  j  redt  +  k0  (0) , 

kx  =  a2j  re(X)dt  +  ^(0) , 

q0  =ri\r®rdt  +  q0(0),  (4.18) 

ql=y2\r®r'dt  +  ql($), 
q2  =r3\r®r"dt  +  q2(0), 

where  Sir  ah  y,  -  positive  adaptation  gains  selected  by  the  system  designer. 

Results  of  the  above  mathematical  analysis  make  it  evident  that  this  approach  does  not 
require  knowledge  of  Omni-Wrist  III  dynamics.  In  addition,  there  is  no  explicit  definition  of  a 
reference  model  to  specify  the  desired  behavior  of  the  system.  However,  signal  0r  applied  to  the 
input  represents  the  desired  dynamics  of  the  system  response;  hence,  this  signal  could  be 
generated  by  a  reference  model  Gm  selected  to  satisfy  the  design  specifications. 

4.2.2.  System  Implementation.  Application  of  the  method  of  Lyapunov  functions  results 
in  a  highly  robust  controller  design.  However,  before  proceeding  with  a  prototype 
implementation  a  couple  of  important  issues,  pertaining  to  system  stability  and  performance, 
should  be  discussed.  Note  that  the  control  law  defined  by  (4.6)  generates  signal  T,  while  the 
physical  input  to  the  dynamic  channel  is  Uj.  A  correspondence  between  the  two  signals  is 
established  by  (4.4),  hence 


17 


Ui=N7\s)Ti  (4.19) 

It  appears  that  by  applying  T,  rather  than  Ut  the  residual  signal  (Nu (s)-l)Tj  is  ignored 
[11].  Even  if  this  signal  is  regarded  as  a  part  of  component  d{,  it  cannot  be  considered  “slowly 
time-varying,”  since  it  includes  signal  7)  that  has  high  frequency  content.  On  the  other  hand,  A)  is 
a  constant  coefficient,  and  it  could  be  demonstrated  that  even  if  equation  (4.5)  is  scaled  by  a 
factor  of  Ni  and  the  proposed  controller  is  described  by  (4.6),  we  could  still  obtain  the  same 
expressions  for  the  conventional  controller  (4.17)  and  adjustable  gains  (4.18). 

Another  problem  could  be  encountered  because  the  cross-coupling  effects  from  the 
actuator  inputs  to  the  azimuth/declination  outputs,  included  in  component  d,  are  very  strong. 
Indeed,  changing  only  one  output  coordinate  requires  the  motion  of  both  linear  actuators,  and 
therefore,  the  application  of  voltage  signals  to  both  motors.  The  adaptive  algorithm  presented  in 
this  paper  is  adequately  suited  only  for  loosely  coupled  systems  [12],  [13];  therefore,  additional 
steps  must  be  taken  to  reduce  the  coupling  effects.  A  decoupling  fdter  must  be  introduced  in  the 
input  of  the  system.  This  filter  is  based  on  the  solution  of  the  inverse  pose  kinematics  problem 
[10]  and  transforms  desired  azimuth/declination  angles  into  corresponding  linear  actuator 
coordinates.  Configuration  of  the  entire  decentralized  system  is  presented  in  Fig.  4.4. 


AFF  -  adaptive  feedforward,  AL  -  adaptation  law. 

The  linear  actuators,  represented  in  the  above  figure  as  DYNAMICS,  take  voltage  as 
inputs  and  provide  actuator  position  as  outputs,  generated  by  optical  incremental  encoders,  which 
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is  the  reason  why  the  controller  are  built  around  the  actuators.  The  KINEMATICS  block 
represents  the  kinematical  structure  of  the  device  coupling  the  outputs,  while  the  INVERSE 
KINEMATICS  FILTER  transforms  the  reference  azimuth  and  declination  coordinates  into 
reference  actuator  positions.  Each  channel  is  controlled  individually,  thus  facilitating 
decentralized  operation. 

4.3.  Bragg  Cell  Control  System. 

A  typical  configuration  of  an  acousto-optic  tracking  system,  such  as  the  one  shown  in 
Fig.  3.2,  includes  two  Bragg  cells  required  to  perform  2-dimensional  beam  steering  and  a 
quadrant  detector,  which  provides  beam  position  feedback  to  the  controller  that  regulates 
frequencies  of  the  RF  signals. 

The  dynamics  of  a  Bragg  cell  is  characterized  by  a  first-order  transfer  function  given  by 
(3.5),  which  is  also  supported  by  the  results  of  our  step  response  experiments  [6],  [7], 
Considering  that  the  access  time  of  these  devices  could  easily  be  on  the  order  of  tens  of 
microseconds  or  less,  their  steering  bandwidth  is  typically  very  large  (usually  on  the  order  of  tens 
of  kHz).  Therefore,  a  simple  gain  controller  in  the  feedback  appears  to  be  sufficient  to  reject 
most  of  the  distortions,  and  an  equation  for  the  control  effort  applied  to  a  Bragg  cell  could  be 
written  as  follows 

f=fc+H*vaz,el,  (4.20) 

where  vaZyei  -  azimuth  or  elevation  feedback  signal  from  the  quadrant  detector. 

This  approach,  however,  does  not  work  in  practice.  Any  quadrant  photodiode  will  act  as  a 
source  of  several  types  of  noise,  including  signal  shot  noise,  background  noise,  and  dark  current 
noise;  while  thennal  noise  will  be  generated  in  the  electronic  circuitry.  System  perfonnance  will 
be  affected  by  all  noise  frequencies  within  the  passband  of  the  tracking  system,  which  will  pose  a 
significant  problem.  Indeed,  a  device  as  agile  as  a  Bragg  cell  would  respond  to  almost  any  signal 
coming  from  a  quadrant  detector,  regardless  of  whether  the  signal  represents  an  actual 
displacement  of  the  laser  beam  or  just  the  additive  noise.  Therefore,  a  constant  gain  controller  in 
the  feedback  needs  to  be  complemented  by  intelligent  filtering  of  the  position  measurement 
signal.  A  block  diagram  of  the  proposed  control  system,  per  channel,  either  azimuth  or  elevation, 
is  presented  in  Fig.  4.5.  In  this  example  it  is  assumed  that  center  frequency  of  the  Bragg  cells  is 
24MHz. 
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Figure  4.5.  Control  system  configuration 

A  disturbance  signal  with  a  specific  spectrum,  e.g.  representing  aircraft  jitter, 
continuously  affects  the  pointing  direction  of  our  transmitter.  If  this  disturbance  is  not 
completely  compensated  by  a  fast  acousto-optic  steering  (FAOS)  device,  the  resultant  pointing 
angle  error  causes  response  in  the  quadrant  detector.  Since  a  signal  from  the  detector  is 
contaminated  with  noise,  it  is  first  filtered,  and  then  used  by  a  constant-gain  controller  to  adjust 
the  frequency  of  the  FAOS  around  /j=24MHz.  The  purpose  of  a  Kalman  filter  is  to  estimate  the 
state  of  a  system  from  measurements,  which  contain  random  errors  due  to  the  noises  in  the 
tracking  system. 

Since  the  controlled  plant  (acousto-optic  device)  in  this  case  is  a  first-order  system,  the 
implementation  of  a  first-order  Kalman  filter  would  probably  be  sufficient  for  rejecting 
measurement  noise.  Generally,  a  first-order  system  could  be  expressed  in  the  discrete-time 
domain  as  follows: 

yn=xn+a-yn- 1.  (4-21) 


where  x  is  the  filter  input,  y  -  filter  output,  and  a  -  parameter  of  the  model. 
Then  an  equation  for  a  first-order  Kalman  filter  is 


^ n\n-\  (  \  , 

yn  = - — - \Xn-a-yn-\)+a-yn_x, 

s  +  Mn  |„_i 


where  s  is  the  noise  variance  and  the  adaptation  mechanism  of  the  filter  is  given  by 


M„ ,  =M„ 

,  •  a 1  +  (x„  -  a 

n\n— 1  n— 

1-M  , 

=  Mn\n  , 

n\n— 1 

n  n\n— 1 

(4.22) 


(4.23) 
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The  above  control  and  adaptive  filtering  algorithms  are  iterative  and  could  be 
programmed  in  software.  The  frequency  at  which  the  control  inputs  to  the  FAOS  system  are 
updated  could  be  very  high  (tens  of  kHz  or  more)  and  is  simply  a  function  of  the  hardware 
characteristics  and  the  latency  in  software-hardware  interaction. 

4.4.  Fusion  of  the  Technologies. 

The  development  of  a  hybrid  beam  steering  system,  combining  the  advantages  of 
different  technologies  fused  together  by  advanced  controls,  is  the  approach  that  has  a  great 
potential.  It  is  well  understood  that  we  need  to  avoid  the  situation  when  both  systems  work  “one 
against  another”  causing  unnecessary  motion,  power  losses  and  resulting  in  positioning  errors. 
The  proposed  architecture  is  presented  in  Fig.  4.6  below.  An  optical  transceiver  with  a  fast 
steerer,  such  as  a  pair  of  Bragg  cells  could  be  installed  on  a  sensor  mount  of  the  Omni-Wrist 
gimbal. 


LOS 


Figure  4.6.  Hybrid  system  configuration 

Coordinated  operation  of  the  proposed  system  could  be  assured  by  hierarchical  control. 
As  could  be  seen  in  the  above  figure,  command  inputs  are  only  applied  to  the  gimbal,  while  the 
line-of-sight  (LOS)  error,  detected  by  a  position  sensing  detector,  is  only  used  to  control  the  fast 
steerer.  This  approach  facilitates  distribution  of  tasks  between  the  two  devices  in  a  very  efficient 
way,  avoiding  the  situation  when  both  are  driven  with  the  same  inputs  and  could  potentially 
work  “one  against  the  other.”  It  is  also  an  efficient  approach  in  the  sense  that  command  signals 
are  typically  generated  to  compensate  for  certain  mechanical  motion  (e.g.  the  changing  attitudes 
of  both  communicating  platforms),  and  are  usually  “slow-speed-large-magnitude”  signals. 
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Contrary,  the  LOS  disturbance,  that  a  position-sensing  detector  is  able  to  pick,  is  a  low- 
magnitude  signal,  which  could  potentially  have  very  high  frequency  content  (if  the  pointing 
distortion  occurs  due  to  high-frequency  jitter  or  the  effect  of  the  optical  turbulence  combined 
with  fast  motion  of  both  the  transmitter  and  the  receiver).  These  intrinsic  properties  lead  to  an 
intuitive  conclusion  that  most  of  the  attitude-related  changes  should  be  handled  by  a  gimbal, 
while  the  fine  LOS  tracking  should  be  primarily  assigned  to  a  fast  steerer.  Additional  interaction 
between  the  two  components  of  a  hybrid  system  is  still  needed,  and  it  will  be  discussed  later  in 
this  section.  Detailed  operation  of  the  system  is  as  follows. 

Command  inputs  in  the  desired  coordinate  system  (yaw/pitch,  azimuth/declination,  etc.) 
are  fed  into  the  Reference  Models  with  desired  dynamics  given  by  the  transfer  functions  Gm, 
producing  the  reference  vector  in  the  desired  coordinate  system.  The  reference  model  vector  0, 
is  transformed  through  the  Inverse  Kinematics  block  (see  Fig.  4.4)  into  the  reference  vector  in 
the  coordinates  of  the  actuators  (typically  counts  of  the  optical  quadrature  encoders).  The 
controller  for  each  axis  converts  the  state  error  (difference  between  the  state  of  the  Reference 
Model  in  the  actuator  coordinates  and  the  state  of  the  plant)  into  a  control  signal,  facilitating  an 
exponential  decay  of  the  state  error  signal  with  predefined  rate,  and  consequently,  decoupled 
operation  in  the  reference  coordinates  [14].  Although  the  resultant  control  law  seems  to  be 
formidable,  a  dedicated  modem  PC  implements  it  quite  effortlessly,  performing  the 
computational  cycle  on  the  20  kHz  clock.  The  main  task  of  the  coarse  steerer  could  be  identified 
as  platform  stabilization  achieved  by  compensating  for  the  disturbances  in  the  orientation  of  the 
optical  platform  as  soon  as  they  are  detected.  Detection  of  such  disturbances  could  be  performed, 
for  example,  by  an  inertial  measurement  unit,  consisting  of  a  3-axis  accelerometer,  a  3-axis 
magnetometer,  and  a  3 -axis  gyroscope,  also  known  as  the  MARG  (Magnetic/ Angular  Rate/Gyro) 
sensor  operating  in  tilt-compensated  (strap-down)  mode.  The  accurate  estimation  of  the 
disturbance  signal,  used  for  the  compensation,  could  be  achieved  by  the  application  of  a  Kalman 
filter  fusing  the  data  collected  from  an  inertial  measurement  unit.  However,  such  work  is  beyond 
the  scope  of  this  project. 

While  the  position  control  system,  outlined  above,  is  necessary  to  compensate  for  the 
disturbance  in  the  orientation  of  the  platform,  optical  tracking  is  required  to  maintain  a  robust 
link  throughout  a  communication  session.  The  main  purpose  of  this  system  is  to  minimize  the 
receive  power  losses  caused  by  the  uncertainty  in  the  line-of-sight  (LOS)  direction  due  to 
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reference  frame  errors,  errors  in  the  pointing  mechanism  (boresight  errors),  dynamic  changes  of 
the  LOS  due  to  transmitter/receiver  motion,  and  beam  wonder  caused  by  the  optical  turbulence. 
Residual  signals  ei  and  02,  representing  pointing  disturbance  to  the  LOS,  must  be  eliminated  by 
the  fast  steerer.  Position  measurements,  performed  by  a  quadrant  detector,  are  always 
contaminated  with  photodetection  and  thermal  noises;  therefore,  these  signals  are  passed  through 
a  Kalman  filter  prior  to  being  applied  to  the  Bragg  cells,  as  shown  in  Fig.  4.5.  The  errors  are  used 
to  adjust  the  ultrasonic  frequencies  of  the  devices  around  fc=24MHz  to  perfonn  agile  beam 
steering. 

It  should  be  well  understood  that  the  driving  mechanisms  for  the  robotic  manipulator 
must  be  controlled  not  only  by  the  platform  stabilization  system,  but  also  by  the  optical  tracking 
system.  Command  Signal  Generator,  shown  in  Fig.  4.6,  is  responsible  for  computation  of  the 
gimbal’s  command  inputs  based  on  both  the  inertial  measurements  and  the  pointing  errors 
measured  by  the  quadrant  detector.  The  first  component  is  necessary  to  assure  that  the  gimbal 
compensates  for  the  changes  in  the  orientation  of  its  own  optical  platform,  while  the  second 
component  is  required  to  track  the  motion  of  the  other  communication  terminal,  which  could 
very  well  be  outside  the  range  covered  by  the  fast  steerer.  Therefore,  the  function  of  the 
Command  Signal  Generator  could  be  expressed  in  the  fonn 


(4.24) 


where  ©n  and  ©a  are  angular  estimates  from  the  inertial  sensors,  ei  and  62  are  the  pointing 
errors,  is  the  threshold  determined  by  a  practical  range  of  the  fast  steerer,  and  function  /  is 
defined  as  follows 


{  x,  if  x>  0 
[0 ,  otherwise 


(4.25) 


The  threshold  errors  in  (4.24)  should  be  chosen  carefully  to  avoid  a  situation  when  the 
Bragg  cells  could  reach  their  saturation  limits  and  not  be  able  to  compensate  for  the  low- 
magnitude-high-frequency  angular  perturbations  to  the  LOS  in  both  azimuth  and  elevation 
channels. 

4.5.  Simulation  Results. 

The  designed  hybrid  steering  system  has  been  tested  under  various  operating  conditions, 
which  could  arise  from  motion  of  the  communicating  platforms,  as  well  as  the  adverse  effects  of 
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the  mechanical  vibrations  (jitter)  or  the  effects  of  the  optical  turbulence,  resulting  in  the  change 
of  the  beam  direction.  The  results  of  this  simulation  study,  conducted  using 
MATLAB/SIMULINK  models  of  the  system  components,  are  summarized  below. 

First,  the  gimbal  control  system  was  designed  and  tested.  The  selected  goal  was  to 
achieve  the  settling  time  Tset  =  50ms  with  no  overshoot,  which  places  the  poles  of  the  closed-loop 
system  at  -80±  j,  resulting  in  the  following  transfer  function  of  the  reference  model  (see  Fig.  4.4) 


r  (  t  6401 
",{S  s2  +1805  +  6401 


(4.26) 


Fig.  4.7  features  response  of  the  gimbal  to  a  square-wave-shaped  command  applied  to  the 
azimuth  input.  The  adaptation  effects  are  clearly  seen  during  each  excursion  of  the  output 
signals,  which  track  the  desired  trajectory  closer  as  time  goes  on. 
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Figure  4.7.  Response  of  the  gimbal  control  system  to  a  square  wave  signal  applied  to  the  azimuth  channel 
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Figure  4.8.  Response  of  the  gimbal  control  system  to  a  square  wave  signal  applied  to  the  elevation  channel 


Similar  results  could  be  observed  when  the  elevation  command  is  changed  as  shown  in  Fig.  4.8. 

Next,  a  hybrid  steering  system  was  assembled  in  the  simulation  environment  according  to 
the  architecture  in  Fig.  4.6.  In  addition  to  the  control  commands,  used  in  the  previous  two 
simulations,  the  outputs  were  contaminated  with  zero-mean  additive  noise,  emulating  the  effects 
of  the  pointing  errors  caused  by  the  platform  jitter  and  the  optical  turbulence.  To  represent  a 
realistic  environment,  perturbations  to  the  LOS  in  both  azimuth  and  elevation  channels  were 
chosen  to  have  Gaussian  PDFs,  thus  resulting  in  the  radial  pointing  errors  having  Rayleigh 
distribution.  An  important  aspect  of  a  tracking  experiment  is  the  choice  of  the  spectral 
characteristics  of  a  disturbance  signal.  It  could  be,  for  example,  recorded  vibration  spectra  from 
satellites,  aircraft,  ground  vehicles,  etc.  All  these  characteristics  have  very  particular  shapes, 
possibly  with  resonant  peaks,  representing  operation  of  specific  subsystems  onboard  the 
communication  platform  as  well  as  its  motion  patterns. 

For  the  purpose  of  our  study,  we  chose  a  more  generalized  spectrum.  A  disturbance 
signal  was  formed  by  filtering  random  noise  with  a  second-order  low-pass  filter  with  a 
bandwidth  of  2kHz.  This  results  in  almost  flat  spectrum  extending  to  2kHz,  which  exceeds  the 
effects  of  most  of  the  realistic  environments,  where  precise  pointing  of  a  laser  beam  is  adversely 
affected  by  vibrations  and  atmospheric  effects.  A  sample  of  the  temporal  data,  representing  the 
jitter  affecting  the  optical  platform,  is  shown  in  Fig.  4.9.  Also  shown  in  the  figure  below  is 
response  of  the  hybrid  steering  system  compensating  the  effects  of  the  jitter. 


Figure  4.9.  Temporal  response  of  the  system  to  high-frequency  jitter  without  compensation  and  with  hybrid  tracking 
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Fig.  4.10  features  the  spectral  content  of  the  jitter,  shown  in  a  solid  line,  as  well  as  the 
frequency  characteristics  of  the  designed  system  compensating  for  the  effects  of  the  jitter  (dash- 
dotted  line). 


Finally,  system  performance  was  tested  when  control  commands  are  applied  to  the  inputs,  while 
the  outputs  are  subjected  to  the  jitter  noise.  These  results  are  presented  in  Fig.  4.1 1  and  Fig.  4.12 
for  azimuth  command  and  elevation  command,  respectively. 


Figure  4. 1 1.  Response  of  the  hybrid  control  system  to  a  square  wave  signal  applied  to  the  azimuth  channel 
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Figure  4.12.  Response  of  the  hybrid  control  system  to  a  square  wave  signal  applied  to  the  elevation  channel 


It  is  evident  from  the  above  figures  that  the  adaptation  process  is  accelerated  and  the  tracking 


capabilities  of  the  hybrid  system  are  noticeably  enhanced. 

Additional  analysis  of  the  temporal  and  spectral  data  presented  in  Fig.  4.9  and  Fig.  4.10 

3  2 

reveals  that  the  variance  of  the  uncompensated  jitter  is  0.1232*10"  rad“.  When  the  hybrid  system 

5  2 

is  enabled,  the  variance  reduces  to  0.1132*10"  rad“,  or  approximately  by  a  factor  of  108.8.  The 


associated  reduction  of  the  pointing  error  in  the  link  budget  equation  may  be  used  in  several 
ways,  including  the  increase  of  the  link  margin,  extending  the  link  range,  reduction  of  the 
transmit  power,  etc. 
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5.  Additional  Considerations  for  Quantum  Communication  Systems 

Quantum  communication,  a  preferred  modality  when  security  of  a  data  link  is  extremely 
important,  presents  additional  challenges  in  the  development  of  a  pointing,  acquisition  and 
tracking  (PAT)  system.  These  challenges,  not  typically  encountered  in  conventional  laser 
communications,  may  include,  but  not  be  limited  to,  wavelength  issues,  using  separate  sources 
for  tracking  and  communication,  polarization  of  the  transmitted  signal,  etc. 

5. 1.  Wa  ve length  Comp  a  tibility. 

When  very  low  transmit  signal  levels  are  used  to  send  data,  tracking  on  the 
communication  beam  becomes  impossible  and  the  second  laser  source  is  needed  to  send  a 
beacon  signal.  The  extreme  situation  occurs  when  single  photons  are  used  to  encode  the  bits  of 
information.  The  second  laser  source  must  have  a  different  wavelength  from  that  of  the 
communication  signal,  so  that  they  could  be  spectrally  separated  by  the  receiver.  While  sending 
two  aligned  beams  at  different  wavelengths  is  generally  not  a  problem,  in  our  proposed  approach 
it  becomes  a  significant  challenge.  The  Bragg  cells  deflect  each  beam  at  an  angle,  which  is  a 
function  of  the  wavelength,  as  could  be  seen  from  (3.5). 

The  situation  could  be  alleviated  if  the  two  sources  are  chosen  in  such  a  way  that  their 
wavelength  ratio  is  2:1.  As  a  practical  example,  it  should  be  possible  to  use  a  775nm  laser 
source,  and  a  beam  splitter  to  divide  the  source  power  between  two  optical  trains.  Then  in  one  of 
the  trains  we  can  down-convert  775nm  into  1550nm  wavelength,  which  can  be  used  as  a  beacon 
signal  without  having  to  worry  about  its  high  power  levels,  since  this  wavelength  is  considered  to 
be  “eye-safe.”  The  other  wavelength,  775nm,  which  could  be  used  to  encode  binary  information, 
is  not  eye-safe;  however,  the  transmit  power  used  in  quantum  communications  is  always  within 
the  harmless  range. 

Both  wavelengths  could  be  aligned  and  steered  in  exactly  the  same  direction  by  the  Bragg 
cells  if  we  use  the  first  diffracted  order  of  1550nm  to  perform  tracking,  and  the  second  diffracted 
order  of  the  775nm  wavelength  for  communication,  for  which  all  angles  in  (3.5)  will  be  doubled. 
The  only  concern  about  using  the  higher  diffracted  order  for  the  communication  carrier  is  the 
efficiency  and  the  available  power  in  the  transmitted  signal.  To  describe  energy  distribution  over 
the  scattered  orders  of  the  Bragg  cell,  a  set  of  coupled  differential  equations  [15],  [16]  can  be 
considered 
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A 

where  n  is  the  order  of  the  scattered  beam,  a  -  phase  peak  delay  in  the  medium,  andg  is  the 
nonnalized  position  of  the  beam  inside  the  cell. 

To  emphasize  on  the  beam  intensities  on  the  output  of  the  AOD  we  set  E,=  1.  The  phase 


peak  delay  is  defined  as 

A  CkJL 
a  =  — - — , 


(5.2) 


where  km  is  the  light  wave  number,  S  -  acoustic  field  amplitude,  L  is  interaction  length  of  the  cell 
defined  by  the  transducer  size,  and  constant  parameter  C  is  defined  in  [17]. 

When  Q»27t,  the  cell  is  considered  to  be  in  the  Bragg  regime.  In  reality,  Q  is  finite  and 
in  many  instances  can  be  slightly  greater  than  2ti  (this  especially  applies  to  the  case  when  a  large 
steering  range  needs  to  be  achieved).  The  general  trend  is  to  concentrate  the  optical  energy  in  the 
immediate  spatial  neighborhood  of  the  zeroth  order.  Therefore,  truncation  of  much  weaker 
higher  diffracted  orders  creates  negligible  numerical  errors.  Applying  (5.1)  to  the  first  five 
scattered  orders,  we  obtain  the  following  equations  for E_2 ,E_VE0,EVE2  [15]: 


The  set  of  coupled  differential  equations  given  by  (5.1)  can  be  solved  numerically  for  a 
given  number  of  diffracted  orders  to  demonstrate  how  the  change  of  Q  affects  intensity 
distribution  among  them.  The  following  two  figures  show  how  intensities  vary  when  a  Klein- 
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Cook  parameter  is  changed,  and  the  phase  peak  delay  is  assumed  to  be  equal  to  4.  For  the  sake  of 
simplicity  only  six  diffracted  orders  are  used  along  with  the  non-diffracted  order. 

Fig.  5.1  presents  numerical  results  for  Q=2ti.  It  is  seen  from  the  figure  that  there  is  a 
noticeable  amount  of  intensity  in  scattered  orders.  At  the  output  of  the  Bragg  cell  intensity  of  the 
first  diffracted  order  is  approximately  84%  of  that  of  the  incident  light  (which  corresponds  to 
84%  diffraction  efficiency).  The  remaining  16%  is  distributed  among  the  other  diffracted  beams. 
The  situation  presented  in  Fig.  5.1  illustrates  that  by  selecting  a  proper  value  fora  we  can  put  a 
significant  portion  of  the  beam  intensity  into  the  first  order. 


Figure  5.1.  Intensity  distribution  for  Q=2ji 

Fig.  5.2  presents  numerical  solution  results  when  the  Klein-Cook  parameter  is  changed  to 
4n  and  the  systems  is  “deeper”  in  the  Bragg  regime.  As  can  be  seen  from  this  numerical  solution, 
higher-order  diffracted  beams  start  “fading”  compared  to  those  in  Fig.  5.1.  At  the  same  time 
intensity  of  the  first  order  diffracted  beam  at  the  output  of  the  Bragg  cell  practically  does  not 
change  and  stays  around  84%  of  the  incident  light  intensity. 
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Figure  5.2.  Intensity  distribution  for  Q  =  4tc 

Analysis  of  the  above  figures  as  well  as  the  results  presented  in  [18],  [19]  leads  to  a 
conclusion  that  it  is  possible  to  design  an  acousto-optic  beam  steerer  that  will  offer  high 
efficiency  for  the  first  diffracted  order  at  one  wavelength  (1550nm)  and  acceptable  efficiency  for 
the  second  diffracted  order  at  the  other  wavelength  (775nm).  For  single -photon  transmission  the 
output  power  for  the  communication  signal  only  needs  to  be 

Pt  =  hv* B  ,  (5.4) 

where  hv  is  the  energy  of  a  photon  and  B  is  the  bit  rate. 

For  1  Gbit/s  transmission  the  above  value  needs  to  be  as  little  as  2.56*10‘10W  or  -65.9  dBm. 

5. 2.  POLARIZA  TION  COMPA  TIBILITY. 

When  two  orthogonal  polarizations  are  used  as  the  base  quantum  states  for  exchanging 
secret  information,  maintaining  their  orientation  on  the  transmitter  side  is  a  critical  task.  One 
challenge  is  presented  by  the  acousto-optic  Bragg  cells,  which  require  specific  polarization  of  the 
incident  beam  and  rotate  its  plane  by  90  deg.  upon  diffraction,  as  illustrated  in  Fig.  5.3(a).  This 
makes  polarization-based  data  encoding  difficult,  since  this  process  must  precede  the  beam 
steering.  Furthermore,  analysis  of  the  robotic  manipulator,  used  in  our  system,  reveals  that  in  the 
process  of  operation  its  optical  mount  continuously  changes  the  roll  orientation,  as  could  be  seen 
in  Fig.  5.3(b).  Finally,  the  changing  attitude  of  the  communication  platform,  as  a  result  of  vehicle 
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or  aircraft  motion  and  maneuvering,  could  also  rotate  the  plane  of  polarization  of  the  transmitted 
signal,  which  inevitably  leads  to  a  change  in  the  polarization  base. 


Figure  5.3.  Challenges  for  maintaining  orientation  of  the  polarization  state  in  transmitted  signals 

(a)  Specific  polarization  of  input  and  output  signals  of  acousto-optic  deflectors 

(b)  Varying  roll  angle  of  gimbal’s  pointing  mechanism  (cross-shaped  optical  mount) 


The  first  problem,  associated  with  the  use  of  the  acousto-optic  devices,  could  be  solved 
by  implementing  a  polarization  diversity  steerer  featured  in  Fig.  5.4.  This  is  a  relatively  simple 
modification  of  the  optical  setup,  which  requires  twice  as  many  acousto-optic  devices,  but  works 
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Figure  5.4.  Acousto-optic  system  utilizing  diversity 
approach  to  steer  a  beam  with  arbitrary  polarization 


for  an  arbitrary  polarization  of  the 
incident  signal.  To  address  the  second 
challenge  outlined  in  this  section  and 
maintain  integrity  of  the  base  quantum 
states  it  may  be  necessary  to  develop  a 
technology  for  polarization  tracking.  To 


implement  such  system,  mathematical  modeling  techniques  and  sensory  data  of  the  PAT  system 


must  be  used  for  calculating  the  control  input  signal  of  the  polarization  rotator. 


The  proposed  hybrid  steering  device  relies  on  the  information  from  the  MARG  sensor, 


mentioned  in  Section  4.4,  which  provides  complete  information  on  the  orientation  of  the 


communication  platform.  This  information  could  also  be  used  to  compute  rotation  of  the  linear 
polarization  planes  associated  with  the  changing  attitude  of  the  ground-based  or  the  airborne 
vehicle.  Furthemore,  additional  rotation  induced  in  the  beam  steering  process  performed  by  the 
gimbal  could  be  found  by  analyzing  the  kinematics  of  the  device  and  developing  a  mathematical 


model  relating  gimbal  control  inputs  to  the  resulting  pitch,  yaw,  and  roll  coordinates  of  the 
optical  mount.  The  concept  of  the  proposed  system  is  illustrated  in  Fig.  5.5.  Precise  orientation 
of  the  polarization  planes  of  the  transmitted  signal  could  be  maintained  in  real  time,  considering 
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Figure  5.5.  System  for  real-time  compensation  of  polarization  base  distortions 
that  currently  available  electro-optical  rotators  offer  adequate  dynamic  perfonnance  with  respect 

to  the  time  scale  of  mechanical  motion. 

While  most  of  the  approaches  to  make  the  proposed  hybrid  technology  compatible  with 
quantum  communication,  as  discussed  in  this  chapter,  constitute  possible  future  research;  the 
development  of  a  polarization  tracking  system  is  within  the  framework  of  this  extension  project. 
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6.  Polarization  Control 

Polarization  control  system  has  two  major  components:  one  associated  with  the 
estimation  of  the  communication  platform  attitude,  which  could  change  as  a  result  of  motion,  and 
another  one  addressing  rotation  of  the  moving  platform,  which  performs  pointing  of  the  optical 
telescope. 

6.1.  Platform  Attitude  Estimation. 

6.1.1.  Inertial  Sensors.  The  PAT  system  is  required  to  compensate  for  the  vibrations 
applied  to  the  optical  platform  while  the  air  or  ground  vehicle  is  in  motion.  The  degradation  of 
the  performance  of  the  communications  system  is  mitigated  through  the  proper  application  of 
advanced  control  laws.  For  example,  additional  feed-forward  vibration  rejection  control  system 
[20]  could  be  used  that  utilizes  a  set  of  inertial  navigation  sensors  to  measure  the  optical  platform 
orientation  disturbance  and  calculates  the  control  effort  that  drives  the  actuators  of  the  Omni- 
Wrist  III.  The  signals  from  the  inertial  navigational  unit,  consisting  of  a  3-axis  gyroscope,  3-axis 
accelerometer,  and  a  3-axis  magnetic  sensor,  are  ‘fused’  to  form  a  quaternion  representation  of 
the  orientation  of  the  optical  platform.  The  control  system  compensates  for  the  disturbance  in  the 
orientation  of  the  optical  platform  as  soon  as  it  is  detected  due  to  the  feedforward  mode  of 
operation  (in  a  feedback  configuration,  such  as  the  optical  tracking,  the  disturbance  is  rejected 
after  it  degrades  the  performance  of  the  communications  system).  The  development  of  an 
Extended  Kalman  filter  ‘fusing’  the  inertial  navigation  sensor  data  is  presented  below. 

The  performance  of  any  vibration  rejection  system,  and  consequently  the 
communications  system,  relates  directly  to  the  disturbance  measurement.  Within  the  scope  of 
another  project,  an  inertial  measurement  unit  was  designed  utilizing  a  3-axis  accelerometer,  a  3- 
axis  magnetometer,  and  a  3-axis  gyroscope.  The  MARG  (Magnetic/Angular  Rate/Gyro)  sensor 
operates  in  tilt-compensated  (strap-down)  mode. 

The  magnetic  measurements  are  collected  from  a  Honeywell  HMC1023  3-axis  magneto¬ 
resistive  sensor.  The  acceleration  measurements  are  facilitated  by  a  two-axis  Analog  Devices 
ADXL203  MEMS  sensor  and  a  perpendicularly  mounted  Analog  Devices  ADXL103  single-axis 
MEMS  sensor.  The  angular  rate  information  is  provided  by  three  perpendicularly  mounted 
Analog  Devices  ADXRS150  MEMS  sensors.  The  use  of  three  types  of  sensors  is  justified  by  the 
unacceptable  errors  limiting  the  usability  of  an  inertial  navigation  system  equipped  with  a  single 
type  or  any  combination  of  two  types  of  sensors.  The  magnetic  sensors  in  conjunction  with  the 
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accelerometers  provide  orientation  information  that  is  flawed  in  the  presence  of  acceleration 
other  than  gravity  and/or  in  the  presence  of  a  magnetic  field  other  that  of  the  Earth,  while 
orientation  integrated  from  the  angular  rates  measured  by  the  gyroscopes  is  exposed  to  errors 
originating  form  the  zero  drift  of  the  sensors  [21].  The  scenario,  in  which  the  errors  originating 
from  different  sources  distort  the  orientation  vector  provided  by  the  accelerometers  and  magnetic 
sensors  and  the  orientation  vector  calculated  from  the  data  collected  from  the  gyroscopes,  calls 
for  the  application  of  the  Kalman  filter  and  ‘fusing’  of  data  collected  from  all  the  sensors  into  a 
single  orientation  vector. 

6.1.2.  Quaternions.  Reliable  operation  of  the  orientation  measurement  subsystem  is 
dependent  on  proper  representation  of  the  state  of  the  system.  Quaternion  representation,  widely 
used  in  navigation  systems  and  computer  graphics,  was  selected  due  to  its  singularity-free 
characterization  of  orientation.  Basic  concepts  related  to  manipulation  with  quaternions  utilized 
in  this  project  are  summarized  below: 

Quaternions,  numbers  with  three  imaginary  parts,  can  be  described  as  an  extension  of 
complex  numbers  (with  only  one  imaginary  part).  Similarly  to  a  complex  number  c  =  a  +  ib 
given  by  two  real  numbers  a  and  b,  where  i  is  the  imaginary  unit  defined  as  i  =  -1,  a  quaternion 
q  =  qw  +  iqx  +  j  qv  +  k q:  is  given  by  four  real  numbers  qw,  qx,  qy,  and  qz  and  three  imaginary  units 
defined  asi  =  j  =  k  =  ijk  =  -1 .  Quaternions  are  well  suited  for  representation  of  spatial  rotation 
thanks  to  their  compact  notation,  relatively  simple  operators,  and  the  singularity-free  property 
[22],  A  rotation  around  an  axis  specified  by  vector  v  defined  as 

v  =  [vx  vv  vJT  (6.1) 


by  angle  0  is  expressed  in  quaternion  form  as 


f  f  Q\  f  Q\  f  Q\ 

q=cos  —  +  i  sin  —  v+jsin  —  v+ksin  —  v_ 

2  2  \2  y  1 2  J 


Two  consecutive  rotations  corresponding  to  quaternions  qi  =  qw\  +  \qx\  +  \qy\  +  kqz]  and 
q2  =  q wi  +  iqx2  +  \qyi  +  kq:o  can  be  represented  by  quaternion  q  =  qw  +  \qx  +  j  qy  +  key-  given  by 
quaterion  product  q  =  q  i  q2  calculated  as 
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(6.3) 


tfw  *Zwl*Zw2  x\tf  x2  Q-yxQyl  Q  :\Q  zl 
<1 lx  =  +  Vx^wl  +  qylQzl  -  V-A^yl 

qy  =  q^qyi  +qyi9w2  +^-1^2  -^1^2 

q*  =  qwlqs2  +  q  zi^w2  +  q^yi  ~qyiqx2- 

In  this  paper,  numerical  integration  of  the  3-axis  angular  velocity  signal  into  quaternion-specified 
orientation  is  implemented  using  quaternion  multiplication.  If  the  current  pose  is  represented  by 
quaternion  q/  and  the  angular  displacement  during  the  last  iteration  is  given  by  quaternion  q (0k 
defined  as 


=  cos(^|  \o)k  \TS  j  +  i  sin(j  |  cok  \TS  j  +  j  sin^  \o)k  |rs  j  ^  +  k  sin^  \ak \TS  j  ^ ,  (6.4) 

where  Ts  is  the  sampling  time  and 

K  |  =  yj  0)lk  +  0)1yi  +  (6-5) 

is  the  magnitude  of  the  angular  rate  calculated  from  the  constituting  axes  measurements,  then  the 
new  pose  can  be  calculated  as 

q*+i  =  q^q*-  (6-6) 


Equation  (6.6)  facilitating  the  quaternion  multiplication  given  by  (6.3)  can  be  expressed  in 
matrix-vector  fonn  as 
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and  linearized  for  small  values  of  angular  displacement  as 
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(6.7) 


(6.8) 


Equation  (6.8)  can  also  be  expressed  in  a  form  convenient  for  the  formation  of  the  state  transition 
matrix  as 
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(6.9) 


q*+i  -  9/  +G kcok, 

where 

-<lx  -<ly  -<h 

qw  <L  ~  Cly 

-  qz  qw  9x 

qy  -qx  qw  \k 

6.1.3.  Kalman  Filter.  The  Kalman  filter  is  an  optimal  estimator  that  can  be  successfully 
utilized  to  combine  data  with  different  distributions  into  a  state  vector  describing  a  process  with 
a  given  variance.  In  the  following  text,  a  Kalman  filter  is  derived  that  provides  an  optimal 
estimate  of  the  state  vector  from  the  data  collected  from  the  MARG  sensor.  Due  to  the  nonlinear 
relationship  between  the  state  and  the  observation,  an  Extended  Kalman  filter  linearizing  the 
state-observation  relationship  around  the  operating  point  needs  to  be  applied.  The  procedure 
leading  to  the  construction  of  the  Kalman  filter  consists  of  the  following  steps:  First,  the  state 
and  observation  vectors  are  defined.  The  state  of  the  system  is  described  by  the  state  vector  x  as 
x  =  [  A  qy  qz  03 x  (o y  co:]T ,  (6.11) 

where  qw,  qx,  qy,  and  q:  constitute  the  quaternion  q  =  qw  +  i qx  +  j qy  +  kq:  representing  the 
orientation  of  the  moving  (sensor)  frame  with  respect  to  the  stationary  (Earth)  frame,  and  cox,  (Oy, 
and  (oz  are  the  angular  velocities  around  axes  x,  y,  and  z  of  the  moving  frame.  The  data  collected 
from  the  nine  sensors  forms  the  observation  vector 

z  =  [  mx  ™y  mz  ax  ay  a:  gx  gy  g,]T,  (6.12) 

where  mx,  mv,  and  mz  are  the  measurements  from  the  three  axes  of  the  magnetic  field  sensor,  ax, 
ay,  and  az  are  the  measurements  from  the  three  axes  of  the  accelerometers,  and  gx,  gy,  and  gz  are 
the  angular  velocity  measurements  around  axes  x,  y,  and  z  obtained  from  the  three  gyroscopes. 
Then  the  state  transition  matrix  F  defines  the  development  of  the  state  of  the  system  in  time 
through  the  difference  equation 

x*=F*_iX*_1+w*,  (6.13) 

where  w*  is  the  normally  distributed  process  noise  vector  with  covariance  matrix  Q  and  k  is  the 
iteration  index.  The  state  transition  matrix  is  defined  as 


(6.14) 
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where  I3  and  I4  are  3x3  and  4x4  identity  matrices  and  Ga  is  defined  in  (6.10).  The  observation 
vector  /k  is  related  to  the  state  vector  \y  through  the  observation  matrix  Ha  as 


hA+v,, 


(6.15) 


where  va  is  the  normally  distributed  observation  noise  with  covariance  matrix  R.  The  observation 
matrix  can  be  derived  from  the  quaternion  qA  representing  the  orientation  of  the  moving  frame 
with  respect  to  the  stationary  frame  by  exploiting  the  fact  that  this  quaternion  transforms  the 
magnetic  field  vector  in  Earth  coordinates  n/  into  body  coordinates  m;''  (measured  by  the 


magnetometers)  as 


mf  =q 


(6.16) 


while  transforming  the  gravity  vector  in  Earth  coordinates  aE  into  body  coordinates  a®  (measured 


by  the  accelerometers)  as 


<  =q*<(ir1> 


(6.17) 


where  vectors  mf ,  n/  ,  af ,  and  aE  are  represented  by  quaternions  as 


=  [0  i  j  k] 


=  .  0  i  j  k]  ,  a* '  =  [  0  i  j  k]  B  ,  aE  =  [  0  i  j  k] 


(6.18) 


Expressing  the  quaternion  transformation  (6.17)  in  equivalent  matrix-vector  form  as 


B  rp  E 

aA  —  T.  a  , 


(6.19) 


where  the  transfonnation  Ta  is  related  to  quaternion  qA  as 

<il  +  <il  ~  h).  -  ql  2 qxqy  -  2ci-^,  2<ia:  +  2qy<iw 

Ta  =  2 qxqy  +  2 q,qw  q2w  -  q]  +  q)  -  q f  2 qyq,  -  2 qxqw  ,  (6.20) 

_  2qxq,~2qyqw  2qyq.  +  2qxqw  ql-ql~q)  +  q] \k 

and  substituting  for  the  gravity  vector  in  Earth  coordinates  in  (6.19)  as 

af=[0  0-l]T  (6.21) 

yields 

0  1  tyxQz  +  2  qyqw 

af  =  T,  0  =  -  2 qyqz  -  2 qxqw  .  (6.22) 

1  2  2  2.2 

L_1J 

Equation  (6.22)  can  be  linearized  around  the  operating  point  given  by  quaternion  qA-  =  qwk  +  iq.xk 
kyk  k</ZA  as 
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corresponding  to  the  orientation  of  the  moving  frame  in  Earth  coordinates  to  the  acceleration  and 
magnetic  field  vectors  in  the  coordinates  of  the  moving  frame.  Combining  the  above  matrices 
with  the  identity  transfonnation  between  the  angular  velocity  signals  in  both  the  observation 
vector  and  the  state  vector  results  in  the  formation  of  the  whole  observation  matrix  as 
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(6.28) 
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0 

I3 


where  I3  is  a  3x3  identity  matrix.  Having  obtained  a  linearized  state  transition  matrix  in  (6.14) 
and  the  observation  matrix  in  (6.28),  we  can  implement  the  Kalman  filter  as  follows.  The 
iteration  begins  with  the  prediction  stage,  in  which  a  new  state  vector  is  formed  by  applying  the 
state  transition  matrix  to  the  previous  estimate  Xfc-i|*-i  as 

x,  =  Ft_|XA._u_i  (6.29) 


and  the  predicted  estimate  covariance  P*|*-i  is  evaluated  as 

PA|A-1  =  ^Ar-l  -^fe-1  +Qk 

In  the  update  stage,  the  innovation  residual  y k  is  found  as 

with  the  corresponding  covariance  Sa  given  by 

$A  =  HA--1PA|A-1HA-1  +  Ra 

leading  to  the  formula  for  Kalman  gain 

Ka^PaiaXA- 

The  state  estimate  is  then  updated  as 

XA|A  =XA|*-l+KAy* 

followed  by  an  update  of  the  estimate  covariance 
Pi|t=(l7-KiHi)P*|t_I, 


(6.30) 

(6.31) 

(6.32) 

(6.33) 

(6.34) 

(6.35) 


where  I7  is  a  7x7  identity  matrix. 

6. 2.  Sensor  Mo  unt  Roll  Angle  Estima  tion. 

Finding  roll  orientation  of  the  sensor  mount  requires  partial  solution  of  the  inverse 
kinematics  problem,  discussed  in  Section  2  of  this  report.  In  this  case  the  azimuth  and 
declination  coordinates,  known  from  the  control  efforts  generated  by  the  gimbal  control  system, 
are  to  be  used  to  find  the  pitch,  yaw  and  roll  coordinates  of  the  transmitting  telescope. 

The  transfonnation  between  the  stationary  frame  of  the  gimbal  and  the  sensor  mount 
frame  through  Leg  A  and  Leg  B  is  found  using  Denavit-Hartenberg  parameters  [2],  These 
transfonnation  given  by  (2.1)  could  be  expanded  as  follows: 
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A  =  Aj  x  A2  x  A3  x  A4  = 
ClC4(C2C3-S2S3Ca)+ 
(S1S3C4+CIS2S4)Sa  + 
5,54Ca 

S,C4(C2C3-S2S3Cj- 

(CAC.-S^S^,- 

cls4ca 

C4{S2C3+C2S3Ca)~ 

C2S4Sa 
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S\{C2S3  +  S2C3Ca)~ 
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S2S3-C2C3Ca 


cls4{c2c3-s2s3ca)+ 

{SlS3S4-ClS2C4)Sa  -  d(Sl (l - C„ )-  CtS2Sa ) 
5,C4Ca 
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(C1S3S4  +  SlS2C4)Sa  -  d(ci (Ca  -l)-StS2Sa) 
CtC4Ca 

S4(S2C3+C2S3Ca)+ 


C2C4Sa 


dC2Sa 


(6.36) 


and 


B  =  B0xB{xB7xB3xB4  = 

CsCg (C6C7  -S6S7Ca )+  -C5Ss (C6C7  -S6S7Ca  )- 
(S5S7Cs  +C5S6Sg)Sa  +  (S5S7Sg  +  C5S6Cs  )Sa  + 


S5SgCa 


S5CgCa 


-Cg{S6C7  +C6S7Ca )+  Ss{S6C7 +c6s7ca)+ 

C6SgSa  C6CgSa 

S5Cs{C6C7 -S6S7Ca)-  ~S5Ss  (C6C7 -s6s7ca)  + 
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C5SgCa  c5cv 
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C5CgCa 
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dM-Ca)-CsS6Sa) 


-dC6Sa 


d{C5{Ca-l)-S5S6Sa) 

1 


,  (6.37) 


where  a  =  45°. 


Then,  knowing  the  azimuth  and  declination  of  the  device,  pitch  and  yaw  coordinates 
could  be  found  as  per  Equations  (2.3)  and  (2.4). 

The  roll  angle  6n  can  be  found  by  substituting  (2.5)  into 

■S.feS,  +  SjCjCJ-S.CA  =  S„SCS,  +C„C,  (6.38) 

and 

64  (S2C-,  +  C2S3Ca )+  C2C4Sa  =  CoCn  (6.39) 


which  correspond  to  elements  (2,2)  and  (3,3)  in  (2.1),  yielding 
-  ^^(l  +  Cj-QCA  =SaS0S„  +CaCn 
and 

SlS2C2{l  +  Ca)+ClC2Sa=C0C„. 

Adding  (6.40)  and  (6.41)  yields 

sas0s„+cac„+c0c„=  0 

which  can  be  transformed  into 


(6.40) 

(6.41) 

(6.42) 
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S'  Ca+C ' 
Cn  saso 


(6.43) 


or 

tan  M0  =  -COS(J’,M)  +  COS°CT,)  (6.44) 

sin  {pitch)  sin(  vuw) 

Equation  (6.44)  is  consistent  with  the  fact  that  Omni-Wrist  III  is  a  two-degree-of-freedom  system 
supporting  rotations  around  the  pitch  and  yaw  axes  only  -  rotation  around  the  roll  axis  is  given 
by  the  pitch  and  yaw  coordinates. 
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7.  Conclusions 

This  report  presents  a  hybrid  beam  steering  system  for  mobile  quantum  communication 
tenninals.  The  proposed  system  comprises  a  novel  robotic  manipulator  Omni-Wrist  III  that 
provides  a  wide  range  of  motion,  and  an  acousto-optic  Bragg  steerer,  which  assures  the  low- 
magnitude  agile  beam  steering  in  the  vicinity  of  any  point  on  the  operational  hemisphere.  A 
hierarchical  control  system  was  synthesized  to  maximize  the  combined  effect  of  these  devices 
and  utilize  their  advantages  to  the  fullest  extent.  As  could  be  seen  from  the  simulation  results,  the 
hybrid  steerer  enjoys  good  robustness,  while  achieving  high  tracking  accuracy  over  an  extended 
field  of  view.  It  was  demonstrated  that  this  system  facilitates  jitter  rejection  exceeding  lOdB  over 
a  range  of  frequencies  spanning  up  to  a  few  kHz. 

A  number  of  issues  specific  to  quantum  communication  using  polarization-based 
encoding,  and  single -photon  transmission  were  identified.  A  viable  solution  to  the  problem  of 
maintaining  polarization  of  the  transmitted  signal  was  developed.  It  allows  to  identify  rotation  of 
the  polarization  plane  relative  to  the  Earth  coordinates  due  to  changing  orientation  of  the  entire 
communication  platform,  as  well  as  the  steering  performed  by  the  coarse  pointing  device 
(gimbal).  If  necessary,  the  issue  of  polarization  change  in  the  process  of  acousto-optic  steering 
could  be  addressed  in  the  near  future. 

It  may  also  be  necessary  to  perfonn  a  more  detailed  simulation  study  of  the  tracking 
system  performance  in  the  presence  of  other  environmental  effects,  such  as  background 
radiation,  photodetection  noises,  etc. 
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